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The direct transition-matrix approach to determination of the electric polarizabihties 
of quantum bound systems developed in my recent work"'^ is applied to study the 
electric multipole polarizabihties of a two-particle bound complex with a central 
interaction between the particles. Expressions for the electric quadrupole and 
octupole polarizabihties of the deuteron are derived and their values in the case 
of the S'-wave separable interaction potential are calculated. 
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1. Introduction 

The electric polarizabilities of hadronic bound complexes are important physical 
quantities that characterize the deformation properties of composite structures 
containing information on the fundamental nuclear force between their constituents. 

The result of the direct determination of the electric dipole polarizability of the 
deuteron by measuring deviations from the Rutherford scattering for d-^°^Pb elastic 
scattering below the Coulomb barrier, 

aE^B) = 0.70 ± 0.05 fm=^ (Ref.2) , 

is near to the result deduced from data for the total photoabsorption cross section 
using the (T_2 sum rule, 

ag(2H) = 0.61 ± 0.04 fm^ (Ref. 3) . 

The calculated values of the deuteron electric dipole polarizability, performed 
both with the simple separable potential that includes the tensor interaction between 
nucleons^ (a£;(^H) = 0.6311 fm^) and with more realistic nucleon-nucleon potentials^ 
(a£;(^H) = 0.6328 fm^) are within the experimental error of the result of Ref. 3, 
both the data^'^ being nearby. 

The anisotropic character of the deuteron electric dipole polarizability caused 
by the presence of the tensor interaction between the neutron and proton has 
been explored both by the method of the traditional nuclear physics^'^ and in 
the framework of the chiral effective field theory^'^. As it was shown in Ref. 4, 
the longitudinal component (with the electric field along the deuteron spin) a^i 
and the transverse component a%i are essentially differed between themselves: 
a^i = 0.669 fm^ and = 0.555 fm^. (The abovementioned electric dipole 
polarizability aEii^H) is an averaged value of the longitudinal and transverse 
components: aEi = la^i + ^c^ei-) 
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Unfortunately, up to now the deuteron electric polarizabilities of the higher 
multipole order have not been experimentally determined yet. 

The influence of the deuteron polarization (distortion) in the processes of 
scattering and disintegration of the deuteron in the Coulomb field of a nucleus 
at low energies was earlier studied in the papers by Abelishvili and Sitenko^ and by 
Clement 

The aim of this work is to study theoretically the deuteron electric multipole 
polarizabilities by applying the direct transition-matrix approach formulated in the 
recent paper ^. In Section 2 the formula for the electric multipole polarizability of 
the two-body bound system composed of one charged and one neutral particles with 
a central interaction is derived. The use of the t-matrix approach in the case of 
the 5- wave separable interaction potential is carried out in Section 3. Section 4 is 
devoted to the results of calculation of the deuteron quadrupole and octupole po- 
larizabilities and discussion of them. 

2. Transition- matrix approach to determination of the multipole polar- 
izabilities of a two-body bound complex 

The electric multipole polarizabilities of a bound complex aEx are the coefficients 
of the asymptotic expansion of the polarization potential of interaction between a 
charged particle and the center of mass of the complex at distances p much more 
than the complex size, 

KoM - -| E ^ . (1) 

where Cq is the charge of the particle 0. In the case of a two-body complex that 
consists of a charged particle 1 and a neutral particle 2 (for instance, the deuteron), 
the electric polarizability of the multipole order 2^ is given by 

aEx = -2 < V'o I Mxg^{-b)Mx | ^/'o > ■ (2) 

Here ipo is the wave function of the bound complex in the state with the energy —b, 
g^{—b) is the "truncated" Green's operator at the energy —6, 

g'^{s)^Qg{e), g{e) ^ {e - h° - v^2)-\ Q = 1 - P, P = V'o >< V'o I , (3) 

g{e) is the total Green's operator, P is the operator projecting on the ground state 
of the complex, h° and Vi2 are the kinetic energy and interaction potential operators 
of the relative motion of the constituent particles 1 and 2 of the complex, M\ is the 
multipole operator 

= e,r^P^{v, • p) = ei (^) r^P,{v ■ p) , (4) 

\mi2/ 

where ei is the charge of the particle 1 of the complex, is the radius-vector of the 
particle i of the complex relative to its center of mass, ri = ^^r, r is the relative 

771.12 

radius- vector between the particles of the complex, r = ri — r2, p is the relative 
radius- vector between the particle and the center of mass of the complex, is 
the mass of the particle i, = mi+ m2, Px{x) is the Legendre polynomial. The 
unit vectors in (4) are marked by the hat. 
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In the momentum space (with r = itiVk) the function < k | Mx \ iJjq > that is 
contained in Eq. (2) can be written as 

(k I M, I ^o) = i'ei (-I^)\,{k)Px{k ■ p) , (5) 
V mio/ 



where 

/I w \ 

\A,,.A 



(6) 



By its definition (3), the "truncated" Green's operator g'^{e) has no singularity at 
e = —b. To show exphcitly cancellations of singularities in Eq. (3), it is worthwhile 
to express the total Green's function in terms of the transition operator, 

gie)^g''ie)+g°ie)tie)g%e). (7) 

Here g°{e) = {e — h°)~^ is the free Green's operator and t{s) is the transition operator 
satisfying the Lippmann-Schwinger equation 

t(e)^v + vg°{e)t{s). (8) 

In the case of a central interaction potential v, the operator t{e) can be 
decomposed into partial components ti{e), each of them being characterized by a 
definite value of the orbital angular momentum of relative motion of particles I. The 
partial wave expansion of the transition matrix is given by 

oo 

< k|tW(£)|k' >= ^(2/ + l)ti{K k'; s)Pi{k ■ k') . (9) 
1=1 

The 5- wave partial component of the transition operator, tQ{e), has pole singularity 
at the point of the energy of the bound complex e — —b and can be represented as 
the sum of the pole tg and smooth to parts. 



to{e) = t^{e)+io{e) , (10) 

3 >< 

e + b 



p, \uq><Uq\ 

^0 (^) = —T — - (11) 



with the vertex function 

I «o) = I ^o) ■ (12) 

Cancelling the pole singularities of the operators g{e) and Pg{s) aX, e — —b in 
Eq. (3) we write the operator g'^{e) as 

g'^{-b)^g°{-b) - g°{-b)\%){^l;,\-\^|;^){^l;,\g%-b) 
+ g\-b)t{-b)g\-b) . 

Here the operator i{—b) denotes the smooth part of the transition operator t, it is 
the sum of the S'-wave component tQ{—b) and all the higher (with I > 1) partial 
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components of the transition operator t{—b) that also have not singularities at this 
point, 

ii-b) = ioi-b) + t^^\-b) , 

oo 

< k\t^''\-b)\k' >^ ^(2/ + 1)^^^, k'; -b)Pi{k ■ k') . (14) 
1=1 

One can readily see that the smooth 5"- wave operator io{e) at the point £ = —b 
is separable, 

io{-b) =1 uo) {uo \, Ri^ b{uo \ {b + h°)-^ \ uo) , (15) 

with the formfactor in the form of the vertex function (12). 

Substituting the expressions for the "truncated" Green's operator g^{—b) (13) 
and for the smooth components of the transition operator t{—b) and to{—b), (14) 
and (15), into Eq. (2) and taking into account the relations 

{ipo I Ma I Vo) = e25xo , 

(16) 

(V-o I 9°i-b)M, I Vo) = (^0 I Mxg°i-b) I iJo)* = -esf 5ao , 
we write the electric multipole polarizability of the two-body complex in the form 
otEx = - 2(V'o I Mxg°{-b)Mx \ V'o) 

(17) 

- 2(^o|MA(7°(-6)tW(-%°(-6)MA|V^o), A = l,2,3,-- - . 

After integrating in Eq. (17) over angular variables with the use of the formula 
(5) for the function M\ipQ and the partial expansion (14) for the operator t^'''\—b), 
we find the final expression for the electric multipole polarizability of the two-body 
complex with a central interaction between the constituent particles: 



_ 2 2 f ^2^^ \ f°° dkk'^ \ ipx{k) 



_ p p dkedk'k''ipl{k)h{k,k';-b)ipx{k') I 
Jo Jo 47r4 (k^ + b)(^ + b) ■ ^ ^ 



Here, the functions (/Pa, which are built out of derivatives of the wave function of the 
ground state of the complex according to Eq. (6), for the dipole (A = 1), quadrupole 
(A = 2) and octupole (A = 3) polarizabilities take the form 



^\{k) 



-%{k) for A = 1 . 

- -Mk) for A = 2 , (19) 

I -C(fc) + Mik) - for A = 3 , 



and the corresponding partial components of the transition matrix satisfy the 
Lippmann-Schwinger equation 

/•oo dk"k"'^ 1 

tx{K k'; -b) = vx(k, k') - / ^^^vxiK k")^, -Mk\ k'- -b) , (20) 

Jo 27r^ 1- b 
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where Vx{k, k') is the partial component of the interaction potential < kl^lk' >, /ii2 
is the reduced mass of the constituents. 

The formula (18) is rigorous, it provides a possibihty to calculate the multipole 
polarizabilities of both Coulomb and nuclear bound complexes with a central 
interaction between constituents. In accordance with Eqs. (18) and (19), to 
determine the 2'*'-pole polarizability of a two-particle complex, it is necessary to know 
only lower-order derivatives of the momentum-space bound-state wave function and 
a single partial component of the transition matrix at the negative energy of the 
bound state. (The P-, D-, and F-wave partial components of the t-matrix are needed 
to determine the dipole, quadrupole, and octupole polarizabilities, respectively.) 

Compared to the approach, which is traditionally applied to calculate the electric 
dipole polarizability of atomic and nuclear two-body systems using the spectral 
expansion of the Green's operator g{—b) in the complete system of eigenfunctions 
of the Hamiltonian h — h° + Vi2, 

a., = 2i:^?y.5l_^J^, (21) 

where the summation is over all the possible discrete and continuous excited states 
n with the the energy e„ and the wave function ipn, matrix approach (18) is much 
simpler and mathematically clear since it rests only on the bound-state wave function 
and the corresponding partial component of the transition matrix at the negative 
energy of the state, both of them being real functions of relative momentum vari- 
ables. With the use of the t-matrix approach, there is no need for calculation of 
more complicated complex continuum wave functions at all. 



3. 5'-wave separable model 

The expression (18) for the electric multipole polarizability of the two-body bound 
complex is obtained assuming that the interaction between constituents is central. 

A further simplification of the formula (18) for the polarizability of the complex 
can be obtained describing the interaction between the charged particle 1 and 
the neutral particle 2 of the complex with the use of the pure 5'-wave separable 
potential^^ 

Vo{k, k') = -uuo{k)uo{k') . (22) 
In this case the wave function of the two-particle complex equals to 

Mk) ^ ^ . (23) 

Using the Yukawa formfactor (Yu) in Eq. (22) 

and the relationship between the parameters u and /3 that takes place if the potential 
(22) forms one bound state of the complex with the binding energy b = ^j^, 
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we write the normalized wave function of the bound complex as 



In this case, the functions ipx{k) (19) formed from derivatives of the wave function 
(26) are given by 



<^i(A;) = 2Co 



k{2e + /32 + K^) 



2\2 ' 



,^ ... _ k'{3k^ + 3{/3' + K^)k^ + {(5' + P'k' + k')} 

^^('^ = ^ (p + ^.)4(,. + ,.)4 ^ ^ ■ 

Note that the considered model (22) does not take account of the interactions 
in the higher orbital states. Also, the sole non-zero S-wave partial component of 
the transition matrix to is not contained in the expression (18) for the multipole 
polarizabilities. Therefore, in this model the second term of the formula (18), in 
which t\ is present only with A > 0, is equal to zero, and the expression (18) is 
reduced to 

Substituting then the functions (27) into the formula (28) and performing 
integration we obtain the following expressions for the dipole, quadrupole and 
octupole polarizabilities of the two-particle complex: 



1 / ms \ 3 miel 3/3^ + IS^^K + Sl/^^K^ + + 48^k^ + 8^^ 

"^'"12 BHB + kW ' 



= (S) T W^^^^^ ^ ''^'^ ^ ''''''' ^ ''''''' 



12 Vmis/ r P\^ + nf 

^ mie\ 

+6A5/3^k'^ + 960p^K^ + 6A0p^K^ + 332P^k'^ + 96pK^ + 12k^} , (30) 

+4900/3^^2 + 10605/3^/t^ + 18270/3^^^ + 26915/3^^^ + 35840/3^^^ (31) 
+26880/3^/^^ + 18000/3^«;^ + 9760/32«;^° + 3616/3^/^^^ + 800/3«;^2 + 8Qk^^} . 



In the case of the zero-range interaction (/3 oo) the expressions (29) - (31) 
simplify to the form 



- 7- 



aEi 



4 KmuJ 

1 / m2 \ ^ rriiel 1 



1 / m2 \ ^ mic^ 1 



,4 ' 



OiE2 



2 VmisJ ^i^ 
45 / m2 \^ mie^ 1 



.6 ' 



(32) 



a£3 



16 Vmi2y 
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Note that the expression for the electric dipole polarizabiUty of the two-body 
bound complex with the separable interaction (29) is in accordance with the 
corresponding formulae obtained by employing the Dalgar no- Lewis perturbation 
technique^^ in the paper by Friar and Fallicros^^ and in the framework of the three- 
body formalism of the effective interaction between a charged particle and a complex 
in our papers^"^^^^. The t-matrix approach with 5*- wave separable interaction has 
been applied to determine the electric dipole polarizabilities of the triton and lambda 
hypertriton as two-cluster systems in the previous our works"*^^'^^. 

4. Numerical results and discussion 

In this section we apply the formalism developed above to the deuteron as a bound 
complex composed of the proton (the particle 1, — e^, mi = nip, /vripe^ — 
28.819893 fm is the proton Bohr radius) and the neutron (the particle 2, m2 = m„) 
with the binding energy h = K^/2yUp„. We calculate the dipole, quadrupole and 
octupolc polarizabilities of the deuteron both in the model of the S-wave separable 
potential (22) with the Yukawa formfactor (24) and in the zero-range limit (/3 oo). 

We find the parameters of the separable potential (22) and (24), i/ and /3, by 
fitting them to the low-energy p — n interaction data: the deuteron binding energy 



and the triplet p — n scattering length ^Qpn or the asymptotic 5*- wave normalization 
of the deuteron wave function As(^ii), 



^apn = 5.424(3) fm (Ref. 18) , AsCR) = 0.8845(8) fm"^/^ (Ref. 19) . (34) 
In the case of the model (22) and (24) we have 



bCR) = 2.224575(9) MeV (Ref. 17) 



(33) 




2{/3 + Kf 
Pk{2^ + k) 





The results of the calculation are shown in Table 1. The existing experimental 
data for the deuteron dipole polarizability are mentioned in Section 1. 
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Table 1. The electric dipole a^i, quadrupole aE2 and octupole polarizabilities of the 
deuteron, aE\{S), calculated with the use of the separable interaction potential (22) (Yu) and the 
zero-range interaction model (ZR) with the parameters fitting to the p — n data (33) - (34) 



p — n interaction 
[parameters] 


ZR [6(2H)] 


Yu [6(2H), 3ap„] 


Yu [6(2H), As^R)] 




0.3776 


0.6259 


0.6292 


ai<;2('H),fm5 


3.5247 


5.9129 


5.9462 


aE3C'R),W 


92.529 


155.38 


156.26 



It is known^^ that for the determination of the electric dipole polarizability of 
the deuteron the adequate description of the asymptotic behaviour of the deuteron 
wave function is of prime importance, since the weakly bound nucleons in the nucleus 
are with a great probability located at distances large compared to the interaction 
range. For the parameters of the potential fitted to both the p — n bound-state data 
(6(^H) and ^^(^H)), we find in Table 1 the following values of the deuteron electric 
polarizabilities 

OiEiCii) = 0.6292 fm^ , q;b2(^H) = 5.9462 fm^ , and assC^) = 156.26 fm^ . (35) 

The polarizabilities (35) are a little larger (by 0.5 % )than the ones obtained with 
the parameters fitted to the data b and ^Upn (when the potential (22) reproduces a 
somewhat lesser value of the asymptotic normalization: As = 0.8820 fm^"*"^^). 

The deuteron dipole polarizability q;^i(^H) is further changed if the tensor p — n 
interaction is taken into account. Thus, for example, with the use of separable 
tensor interaction potentiaP^'^^ (that reproduces As — 0.8843 fm~^/^) the dipole 
polarizability q;ei(^H) increases to 0.6311 fm^ (Ref. 4) approaching the values 
obtained with the realistic potentials q;£!i(^H) = 0.6328(17) fm^ (Ref. 5). As for the 
inclusion of the interaction in the higher partial states, notice that for the dipole 
polarizability of the deuteron q;ei(^H) it is more important the consideration of the 
tensor p — n interaction (in ^Si +^ Di eigenstate) than that of the interaction in 
the P-wave state^^. The tensor p — n interaction also causes the anisotropy effect 
of the deuteron electric polarization — the dependence of the polarizability a^^CH.) 
on the projection of the total angular momentum (the spin) of the deuteron M that 
describes the direction of the spin relative to the electric field vector. The anisotropy 
effect of the electric dipole polarizability has been previously studied in Ref. 4. The 
application of the i-matrix approach to the study of the anisotropy of the induced 
electric quadrupole moment of the deuteron having also its own quadrupole moment 
(Q = 0.2859(3) fm (Rcfs. 23 and 24)) still remains to be carried out. 

The advantage of the t-matrix formalism are clearly evident by the example 
of its appUcation to the determination of the electric dipole polarizability of a 
simple two-particle system with long-range Coulomb interaction — the hydrogen 
atom. In this case, provided that the ratio of the electron mass to the proton 
mass tends to zero, the exact analytical value of the polarizability is known^^"^^: 
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aEi(ii) = , dB = h'^/m^e'^ is the electron Bohr radius. With this, the formula 
(18) may be directly applied to the hydrogen atom. Substituting into the expression 
(18) the explicit analytical expression for the partial (with I — 1) component of the 
Coulomb transition matrix at negative energy —b obtained previously in Ref. 28 we 
find calculating separately the first (analytically) and second (numerically) terms 
that ^ 

a^i(H) = -a| + 0.321066940 = 0.666831341 A^ = 4.5 a| . 

As distinct from the polarizability of the two-particle nuclear system, for which in 
Eq. (18) the first term containing the free propagator is dominant as compared 
to the second one, the polarizability of the Coulomb system (with the long-range 
interaction) is determined by the sum of the two very nearly equal terms — with the 
free propagator and with the P-wave partial component of the transition matrix at 
the negative energy. 

The developed approach to the problem of the determination of the electric 
multipole polarizabilities, which is based on the two-particle transition matrix, is 
found to be rather effective as compared with the traditional method that uses the 
spectral expansion of the Green's operator of the system. The results by Castillejo et 
aP^ testify that with the use of the spectral expansion all the virtual excited bound 
states contribute 81.4% of the magnitude of the electric multipole polarizability of 
the hydrogen atom, the rest is accounted for by the states of the continuum. With 
the application of th matrix approach the necessity of considering a large number of 
virtual states no longer arises at all. In the case of the dipole polarizabihty (A = 1), 
beside the term with the free virtual propagator, the formula (18) contains only one 
term with the P-wave partial transition t-matrix at the negative energy. Moreover, 
the transition matrix at negative energies is a real quantity, as opposed to more 
complicate complex functions of the continuum, which are used in the traditional 
approach. 

Existing advantages of the i-matrix approach motivate necessity of its further 
development with the aim of applying it for two-particle systems with noncentral 
interactions (the tensor n — p interaction) to study the anisotropic properties of 
the dynamic quadrupole polarization of the deuteron. Also, the extension of the 
transition- matrix formalism to investigate the polarizabilities of more intricate three- 
and four-body nuclear systems based on the Faddeev and Faddeev-Yakubovsky in- 
tegral equations is under way. 
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